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Abstract 

We report a high-statistics measurement of the relative branching fraction B(D° — ► 7r + 7r _ 7r )/ 
B(D° -> R-tt+tt ). A 357 fbr 1 data sample collected with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider was used for the analysis. The relative branching fraction B(D° — > 
7r + 7r~7r )/ B(D° — > K~7r + 7r°) is determined with an accuracy comparable to the latest world 
average value. 

PACS numbers: 13.25.Ft, 14.40.Lb 
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INTRODUCTION 

This measurement is the first step towards a high-statistics Dalitz-plot analysis of the 
D° — > 7r + 7r - 7T° decay. The latter could give insight into the controversy on the S-wave 
7r + 7r~ contribution in these decays 0, 0], as well as a sensitive study of the CP violation 
in the neutral D meson system. Knowledge of B(D° — > p-K)/B(D° — > K*K) (also based 
on the D° — > 7r + 7r~7r° Dalitz analysis) could improve our understanding of the apparent 
discrepancy of the measured two-body branching fractions (D° — > KK, tttc) with the 
theoretical expectations The accuracy of the value of B(D° — > 7r + 7r~7r°) as reported in 
PDG04 0| is poor. Using a large data sample of D° decays accumulated with the Belle 
detector, we provide a significantly improved measurement using the D° — > K~TT + ir° decay 
mode for normalization. Since both decay modes involve a neutral pion and the same 
number of charged tracks in the final state, several sources of the systematic uncertainties 
are avoided in a determination of the relative branching fraction. The obtained result can 
then be compared to recent measurements by the CLEO ji| and BABAR 0| collaborations. 
A detailed study of the D° — > tt + 7i^it decay as well as of other D° CP-symmetric final 
states, can be used to further improve statistics for the measurement of the angle 03 (7) of 
the CKM-matrix . 



EXPERIMENT 

The Belle detector is a large- solid- angle magnetic spectrometer located at the KEKB 
e + e~ storage rings, which collide 8.0 GeV electrons with 3.5 GeV positrons and produce 
T(4S) at the energy of 10.58 GeV. Closest to the interaction point is a silicon vertex detector 
(SVD), surrounded by a 50-layer central drift chamber (CDC), an array of aerogel Cherenkov 
counters (ACC), a barrel-like arrangement of time-of-flight (TOF) scintillation counters, and 
an electromagnetic calorimeter (ECL) comprised of Csl (Tl) crystals. These subdetectors are 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return yoke located outside the coil is instrumented to detect K Q L mesons and identify 
muons. The detector is described in detail elsewhere 0, Ef . 



DATA SELECTION 

For this analysis, we used a data sample of 357 fb~ x accumulated with the Belle 
detector. D° candidates are selected from D* — > D°ti s \ ow decays where the charge of the 
TTsiow tags the D° flavour: D* ± — > D° / D°ir^ ow . D* J s originate mainly from continuum. 
Although we do not apply any topological cuts, the yield of D*'s from B- meson decays is 
negligible; such events are rejected by other kinematical cuts such as the strong p cms (-D*) 
requirement. D° mesons are reconstructed from combinations of two oppositely charged 
pions (or a charged pion and kaon in the case of A^7r + 7r°) and one neutral pion. The 
latter is reconstructed from two 7 candidates satisfying the ir° mass requirement given below. 

The following kinematic criteria are applied to the charged track candidates: the distance 
from the nominal interaction point to the point of closest approach of the track is required 
to be within 0.15 cm in the radial direction (dr) and 0.3 cm along the beam direction (dz). 
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We also require the transverse momentum of the track p± > 0.050 GeV/c to suppress beam 
background. Kaons and pions are separated by combining the responses of the ACC and 
the TOF with the dE/dx measurement from the CDC to form a likelihood £(h) where h 
is a pion or a kaon. Charged particles are identified as pions or kaons using the likelihood 
ratio TZpin = £(K)/(£(K) + £(ir)). For charged pion identification, we require TZpin < 0.4. 
This requirement selects pions with an efficiency of 93% and misidentified kaons with an 
efficiency of 9%. For the identification of charged kaons, the requirement is TZpm > 0.6; in 
this case, the efficiency for kaon identification is 86% and the probability to misidentify a 
pion is 4%. 

We impose conditions on the energies of the photons constituting the ir° candidate 
(£ 7 (7T°) > 0.060 GeV), the two-photon invariant mass (0.124 GeV/c 2 < M{^) < 0.148 
GeV/c 2 ) and the 7r°'s momentum in the laboratory frame (piatX 71 " ) > 0.3 GeV/c) to 
suppress combinatorial 7r°'s. 

The mass difference of D* and D° candidates should satisfy the restrictions: 
0.1442 GeV/c 2 < A/f(7r s i ow 7r + 7r-7r°) - M(7r+7r"7r°) < 0.1468 GeV/c 2 and 0.1440 
GeV/c 2 < M(vr s i ow fs:-7r + 7r )-M(^7r + 7r ) < 0.1470 GeV/c 2 . The momentum of the D* 
in the center-of-mass frame of the T(4S) must lie in the range: 3.5 GeV/c < p C ms{D*) < 4.3 
GeV/c. The lower cut is applied to suppress slow fake D*'s reconstructed from combina- 
torial background that originates from B decays. The upper cut restricts p cms (-D*) to the 
region where the Monte Carlo (MC) distribution is in good agreement with the data (see 
Fig. 1, left). To eliminate background from the D° —* K^ir — > (7r + 7r~)7r° decays, the 
following veto on M 2 (n + 7r~) is applied: 0.21 GeV 2 /c 4 < M 2 {it + tt-) < 0.29 GeV 2 /c 4 . 



EFFICIENCY CALCULATION 

To obtain detection efficiencies, about 1.2xl0 6 phase-space distributed MC events have 
been generated for each of the two modes, processed using the GEANT-based detector 
simulation and reconstructed with the same selection criteria as the data. 

To obtain a more realistic D° decay model than a uniform Dalitz-plot distribution, the 
events of both signal MC samples have been weighted using matrix elements based on 
decay models obtained by CLEO ^ in the framework of a 3-resonance model (p + , p~ , p° 
and a nonresonant contribution) for D° — > 7i + n^n° and in the framework of a 7-resonance 
model (p+ K*~, K°* , A" (1430)-, K (1430Y> p(1700)+, A"*~(1680) and a nonresonant 
contribution) for the case of D° -> K-ir+ir [10j (see Fig. 2). 

The obtained yield is normalized to the same MC data but before detector simulation and 
before application of selection criteria. Since for the relative branching fraction measurement 
only the ratio of 7r + 7r~7r° and K~tt + it efficiencies is needed, we only generate events in the 
p cms (D*) range from 3.5 to 4.3 GeV/c where MC is in good agreement with the data (see 
Fig. 1, left). The obtained efficiencies are: e(7r + 7r"7r°) = (13.433 ± 0.077)%, e(^7r + vr°) = 
(13.065 ± 0.074)%. 
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FIG. 1: Left: 

Pcms (-C*) distribution for experimental data (points with error bars), generic MC 
(solid line), charm (originating from e + e~ — > cc) MC including signal events (singly hatched his- 
togram) and background charm MC (doubly hatched histogram). Right: M(7r + 7r~7r°) distribution, 
MC, fitted with a bifurcated hyperbolic Gaussian and a regular Gaussian. 

BACKGROUND DESCRIPTION 

To describe the background shape in the M(tt + 7t~7t ) signal region, a sample of generic 
MC, equivalent to ~250 fb" 1 , has been processed with the same selection criteria as data. 
The contributions of uu, dd, ss, cc and BB backgrounds have been summed up. For the 
cc sample signal events were excluded. Figure 3 shows the individual contributions to the 
M(7r + 7r~7r°) distribution where also a partial data sample, corresponding to the luminosity 
of MC, is shown. The M(7r + 7r _ 7r°) background distribution is fitted with a 2 nd order 
polynomial multiplied by an error function (kaon misidentification region) plus a l si order 
polynomial (combinatoric background) and a small Gaussian peak in the signal region (see 
Fig. 3, right). The Gaussian contribution is mainly due to combinations of a correctly 
reconstructed D° and a random 7r slow candidate. 

Most of the background is from e + e~ — > cc; the BB background is negligible and uu, 
dd, ss backgrounds are linearly distributed in the M(7r + 7r~7r°) signal region. Among the 
cc background sources D* — > D (K7nr )7r is dominant: charged kaons are misidentified as 
pions and M(7rf a k e 7r7r°) is typically shifted downwards by ~ 0.1 GeV/c 2 thus being well 
separated from the signal. 

The background shape in the M(K~7i + ir°) distribution is obtained using a generic MC 
sample. The level of background in this mode is low but still has a nontrivial structure 
(see Fig. 4, right). The distribution has three distinctive peaks: the rightmost one is due 
to misidentified pions from D° — » 7r + 7r _ 7r°, the central one has the same origin as the one 
in the D° — > 7r + 7r~7r° case (random vr s i ow ), the leftmost feature originates mainly from 
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FIG. 2: Dalitz plot distribution for "initial" signal MC (before skimming and detector simulation) 
in the D° — > 7r + 7r _ 7r° (left, 3-resonance model) and D° — ► il'~7r + 7r° (right, 7-resonance model) 
cases. 




FIG. 3: Left: M(7r + 7r _ 7r°) distribution for experimental data (points with error bars), generic 
MC (solid line), charm (originating from e + e~ — ► cc) MC including signal events (singly hatched 
histogram) and background charm MC (doubly hatched histogram). Right: M(7r + 7r~7r°), "charm" 
MC with the signal excluded, fitted to the function described in the text. 

D° — > Knir where n > 3. The distribution is fitted by the sum of a 2 nd order polynomial 
and a 2 nd order polynomial multiplied by an error function (left peak) and two Gaussians. 
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FIG. 4: Left: signal M(Ktt7t°) distribution fitted with 2 bifurcated Gaussians + Gaussian (signal 
peak) and the generic MC shape (background). Right: background shape M(K7rir°), generic MC. 

SIGNAL FIT 

The shape of the signal peak in the experimental 7r + 7r~7r invariant mass distribution 
(see Fig. 5, left) is partially fixed to the MC one [3]: the latter was fitted with a bifurcated 
hyperbolic Gaussian 

Urn 

and a regular Gaussian (see Fig. 1, right). The obtained shape 
with floating normalization and cr's, together with the background shape with its floating 
normalization (i.e. 5 free parameters) is then used as the fit function for the measured 
M(ir + 7T^7r°) distribution. 

To fit the experimental Kirn invariant mass distribution we use a Kirn background 
shape with floating normalization and the sum of two bifurcated and a regular Gaussian 
for the signal peak (see Fig. 4, right). Here, the parameters of the signal peak are free in 
the fit since the level of background is low and statistics are large enough. 

The yields of signal events in each channel, as obtained from the fit, are given in Table |1J 



TABLE I: Number of signal events and efficiencies 



D° decay mode 


N ev 


e, % 


K~7T + TT° 


22803 ± 203 
237520 ± 500 


13.433 ± 0.077 
13.065 ± 0.074 
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FIG. 5: Left: signal M(7r + 7r _ 7r°) distribution, fitted to the signal MC shape (Fig. 1, right) for the 
signal peak and with the generic MC shape (background). Right: misidentification shape, generic 
MC, used for an alternative background fit to estimate the corresponding systematic error. 



SYSTEMATIC UNCERTAINTIES 

Uncertainty on the tracking efficiency for the two charged tracks - 7r + 7r~ or K~n + - 
cancels in the ratio of D° — > 7r + 7r _ 7r° and D° —>■ K~n + 7T° branching fractions. The same 
holds for 7T° and slow pion (from D*) reconstruction efficiencies. A possible difference in 
the efficiency of particle identification selection criteria between MC and data is taken into 
account as a correction to the branching ratio (r^/r^ = 0.962 ± 0.015) and the uncertainty 
of this correction contributes to the systematic uncertainty of the result. The polar angle 
and momentum dependent data/MC corrections are measured independently using a large 
sample of D° — > K~n + decays. The uncertainty in the correction yields a systematic error 
of 1.6%. 

To estimate the systematic error due to the model used to describe the D° — > 7r + 7r~7r° 
decay, we have compared MC efficiencies using the two different models mentioned above: 
the one obtained by CLEO 0] (used in the final result for the branching ratio) and the same 
model but without interference: only the branching fractions of the intermediate resonances 
have been taken into account. Both models consider three intermediate resonances (p + , p~ , 
p° and a nonresonant contribution). 

For the D° — > Knir mode results were obtained by CLEO 0] in the framework of 
a 7-resonance model (p+, K*~, K°* , K (1430)", J? (1430)°, p(1700)+, K*~(1680) and a 
nonresonant contribution) and the 3-resonance model (K*~ir + , K° 7r°, p + vr~, nonresonant 
K~n + 7r°). The difference in the resulting efficiency using two different models is treated as 
a corresponding systematic error. The resulting systematic error due to the decay model 
uncertainty is 1.8%. 
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As mentioned above, we use the generic MC shape in the M(ir + 7r~7T°) signal region as the 
default background description (when calculating the central value of B(D° — > 7r + 7r~7r )). 
To estimate the systematic error related to the description of background, we use an 
alternative function for the background fit. To do so, we extract true Kirir events from 
generic MC reconstructed as 7r + 7r~7r°. The shape obtained (see Fig. 5, right) is fitted 
with a 2 nd order polynomial multiplied by an error function and a I s * order polynomial. 
We use this parameterization for the misidentification background, while another I st order 
polynomial is added to fit the other charm- (other than misidentification), light- and 
6-quark linear contributions. The difference between this and the default fit (0.4%) is 
added in quadrature to the systematic error. 

To determine the signal fit uncertainty, a different D° — > 7r + 7r~7r° signal parameterization, 
composed of two bifurcated Gaussians and one regular Gaussian is used. For the D° — > 
Kirn decay channel we use a bifurcated hyperbolic Gaussian and a regular Gaussian. The 
relative differences to the default fit are found to be 0.3% and 0.5% for the 7r + 7r~7r° and 
K~7T + 7r° mode, respectively, and are added to the total systematic uncertainty. 

Finally, we varied the selection criteria and estimate the systematic error corresponding 
to the possible inadequate background description. The particle identification selection 
was changed to 7^pid(vt) < 0.2, 7^pid(-£0 > 0.8, and the resulting systematic uncertainty 
is found to be negligible. Change of the p cm s(D*) requirement (3.0 GeV/c) results in a 
0.4% error. A significant uncertainty is due to the Ks veto; after excluding the region 
0.195 GeV 2 /c 4 < M 2 (tt + < 0.305 GeV 2 /c 4 , the resulting tt+tt^tt yield changes by 1.6%. 

The individual sources of the systematic error are listed in Table ITT1 The total uncertainty 
is obtained by adding all contributions in quadrature. 



TABLE II: Systematic errors 



Source 


Error, % 


MC statistics 


0.8 


PID efficiency of K/ir 


1.6 


Decay model 


1.8 


Fit (background & signal) 


0.7 


Pcms(D*) cut 


0.4 


K° s veto 


1.6 


Total 


3.1 



RESULTS 

Summarizing the discussion above, we obtain the following ratio of the branching frac- 
tions: 
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n n TV 



N 



B(D° — > i^~7T + 7T ) N^-^O e^+Tr-TrO 



• — = 0.0971 ± 0.0009 stat ± 0.0030 

rv 



syst- 



We multiply the obtained value by the 2006 world avearge of B(D° — > K~tv + tv°) = 
(14.1 ± 0.5)% jll[ to calculate the absolute branching fraction for D° —> tv + tv~tv° 
decay (see Table lIHj) . In a recent study by CLEO [5j the relative branching ratio 
B(D° -> 7r + 7r-7r°)/i3(D° -> fsTvr) is measured to be 0.344 ± 0.013. The comparison of the 
corresponding values for the absolute branching ratio B(D° — > 7r + 7r~7r°) shows that both 
results are in good agreement within the errors (see Table UTT]) . 



TABLE III: Values of the absolute branching fraction B(D° — > 7r + 7r~7r ) by Belle and CLEO-c. 





N cv . 


B(D° tv+tv~tv°), lO" 3 


Belle 
CLEO-c 


22803 ± 203 
10834 ± 164 


13.69 ± 0.13 stat ± 0.42 syst ± 0.49 norm 

13.2 ± 0.2 sta t ± 0.5 sys t ± 0.2 norm ± O.lcPcorr. 



Finally, we can compare our measurement of the ratio with the ratio obtained from the 
latest world average values of %^ d q*Z}c-J+) an< ^ B(D G — > 7r _ 7r + 7r°) [ll|. Our result (9.71 ± 
0.31)% is consistent with the world average one (9.29 ± 0.54)% and has better accuracy. 
By choosing the normalization mode D° — > K~tv + tv° we avoid many sources of systematic 
uncertainty including the tv° detection efficiency and uncertainty in the tracking efficiency 
The PID efficiency partially cancels out. 



SUMMARY 

Using 357 fb _1 of data collected with the Belle detector, the first direct measurement of 
the relative branching fraction B(D° — > tv + tv~tv )/ B(D° — > K~tv + tv°) has been performed. 
Our preliminary result 0.0971 ± 0.0031 is compatible with the world average 0.0929 ±0.0054 
and is more precise. The Belle result differs by ~ 2.8 a from the value recently obtained by 
BaBar 0]: 0.1059 ± 0.0006 ± 0.0013. The corresponding value of the absolute branching 
fraction B(D° -> tv+tv'tv ) is (13.69 ± 0.66) xlO" 3 . 
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